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Abstract 

In this paper, we address a classical case of the Calderon (or conductivity) inverse problem in 
dimension two. We aim to recover the location and the shape of a single cavity u! (with boundary 
7) contained in a domain O (with boundary T) from the knowledge of the Dirichlet-to-Neumann 
(DtN) map : f 1—where is harmonic in O \ w, = / and u-^\j = c-t, being 
the constant such that f dnU^ ds = 0. We obtain an explicit formula for the complex coefficients 
am arising in the expression of the Riemann map z 1 —> aiz + oq + that conformally 

maps the exterior of the unit disk onto the exterior of ui. This formula is derived by using two 
ingredients: a new factorization result of the DtN map and the so-called generalized Polia-Szego 
tensors (GPST) of the cavity. As a byproduct of our analysis, we also prove the analytic dependence 
of the coefficients Um with respect to the DtN. Numerical results are provided to illustrate the 
efficiency and simplicity of the method. 


1 Introduction 

Let R be a simply connected open bounded set in with Lipschitz boundary L. Let u be a positive 
function in L°°{0,) and consider the elliptic boundary value problem: 

—V • (aVu) = 0 in 0 (1.1a) 

u = f on r. (1-lb) 

Calderon’s inverse conductivity problem m can be stated as follows: Knowing the Dirichlet-to- 
Neumann (DtN) map A.^ : / 1 —)• dnU^, is it possible to recover the conductivity cr? 

In this work, we focus on the particular case of piecewise conductivity with infinitely high contrast 
(see for instance Friedman and Vogelius Hal who considered this problem in the case of small inclu¬ 
sions). More precisely, we suppose that D contains a cavity w, where uj is an open connected set with 
Lipschitz boundary 7 and such that uJ C D (see Figure [^. We denote by n the unit normal to F U 7 
directed towards the exterior of D \ dJ. 
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Figure 1: The geometry. 


For every / in Ff 2 (F), we denote by c-^) G {Vl\uj) x M the solution to the Dirichlet problem: 
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where c^ is the unique constant such that: 

[ dnU^ 

dcr = 0 . 

( 1 . 2 d) 
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Problem ( |1.2[ ) is well-posed and its solution is the limit of the solution of (1.1) for piecewise constant 
conductivity, when the contrast between the cavity and the background tends to inhnity (see Propo¬ 
sition |A. 1 | of the Appendix for a precise statement of this classical result and for the proof, which is 
given for the sake of completeness). 

Loosely speaking (the exact functional framework will be made precise later on), the inverse 
problem considered throughout this paper is the following: knowing the Dirichlet-to-Neumann (DtN) 
map : f 1 —> dnU^, how to reconstruct the cavity co? 


Remark 1.1 In dimension 2, it is classical to see as the harmonic conjugate function of , i.e. 
the solution to: 


0 

II 

<1 

1 

in Il\uj 

(1.3a) 

II 

on F 

(1.3b) 

0 

II 

on 7 , 

(1.3c) 

where r := is the unit tangent vector to F. The function 

is usually referred to as the stream 


function associated to the potential function . On F, dnU^ = —drV^ and therefore, the knowledge of 
(i.e. the DtN for ) is equivalent to the knowledge of the Neumann-to-Dirichlet map for . 

Classically for inverse problems, the questions of uniqueness, stability and reconstruction have 
been studied in the literature for cavities identification. Regarding uniqueness, it is well-known that 
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one pair {f,dnU-^) of Cauchy data uniquely determines the geometry of the cavity for a Dirichlet 
boundary condition (see Kress [33]) or a Neumann boundary condition (see Alessandrini and Rondi 
m) For Robin type condition, Bacchelli |9| proved that two excitations /i and /2 uniquely determine 
the cavity provided they are linearly independent and one of them is positive. Concerning stability, 
as shown by Mandache [36] . logarithmic stability is best possible (see also Alessandrini and Rondi 
|2| and references therein). Among the reconstruction methods available in the literature for shape 
identihcation, one can distinguish two classes of approaches: iterative and non iterative methods (see 
for instance the survey paper by Potthast |39j for an overview of reconstruction methods). In the 
first class of methods, one computes a sequence of approximating shapes, generally by solving at each 
step the direct problem and using minimal data (typically only one or several pairs of Cauchy data, 
and not the full DtN map). Among these approaches, we can mention those based on optimization 
unKis], on the reciprocity gap principle [aiiEiiiiii, on the quasi-reversibility mm or on conformal 
mapping [DiMiiiniiiiiiiiiisiiiis]. 


The second class of methods covers non iterative methods which are generally based on the con¬ 
struction (from the measurements) of an indicator function of the inclusion(s). These sampling/probe 
methods do not need to solve the forward problem, but require the knowledge of the full DtN map. 
Among these reconstruction techniques, let us mention -with no claim as to completeness- the enclo¬ 
sure and probe method of Ikehata jSZlESlESlEnidB], Kirsch’s Factorization method dSlElES] and 
Generalized Polya-Szego Tensors in the case of small inclusions igiTllHllS]. 

Our purpose in this paper is to propose a new non iterative reconstruction method that combines 
some of the ingredients used in earlier works, namely: a new factorization result (Theorem 3.1), 
Generalized Polya-Szego Tensors and conformal mapping. The main feature of our reconstruction 
method is that we end up with an explicit reconstruction formula (Theorem |3.4| ) for the complex 
coefficients arising in the expression of the Riemann map z i—)■ aiz -)- uo + that 

conformally maps the exterior of the unit disk onto the exterior of oj. Let us emphasize that these 
reconstruction formulae also yield the analytic dependence of the coefficients with respect to the DtN. 

The proposed reconstruction algorithm can -in principle- be adapted to other boundary conditions. 
However, such as most direct reconstruction methods, it requires the knowledge of the full DtN map 
and so far, it is limited to the two-dimensional case due to the use of conformal mapping. 

The paper is organized as follows: we present in Section]^ a boundary integral formulation of the 
problem. Section]^ is devoted to the derivation of the reconstruction formula, using a new factorization 
result and GPST. Some issues about stability are also discusses therein. Finally, some numerical results 
are given in Section]^ 


2 Boundary integral formulation 

2.1 Background on single layer potential 

In this section, we collect some well known facts of potential theory, and more especially on single 
layer potential, that are crucial for our method. For more details and for the proofs, we refer the 
interested reader to the monographs of McLean m. Steinbach m or Hsiao and Wendland [26] . 
Throughout the article, we shall denote by 

^( 2 :) = -7:^log|x| 

ZTT 

the fundamental solution of the operator —A in M^. 


3 



Let Cj be a bounded, Lipschitz domain (see m Definition 3.28]) and denote by ^ its boundary. 
Let n be the unit normal to ^ directed towards the exterior of Cj. 

The exterior of C* is denoted Ce := \ C*. Given a function u in we denote by Ui and 

Ue its restrictions respectively to Cj and Ce and by [u\<^ = Ui\^ — Ue|<^ G H 2 {'W) its jump across 
We also define similarly the jump of the normal derivative: [dnu]^ = {dnUi)\^ — {dnUe)\^ G 

Definition 2.1 For every q G we denote by the single layer potential associated to the 

density q. 

The single layer potential ^<ifq defines a harmonic function in \ The operator is an integral 
operator with weakly singular kernel, so that for q G for instance and x G \ ^, it reads: 


= / G{x - y)q{y) day. 


Moreover, the single layer potential defines a bounded linear operator from H 2 (“^j into 
and .y<^q admits the following asymptotic behavior at infinity (see for instance ISZl p. 261]) 

y<^q{x) = log |x| + 0(|x|"^), (2.1) 

where (•, •)_! 1 ^ stands for the duality brackets between and This shows in partic¬ 

ular that ^<^(1 G Lf^(M^) if and only if q belongs to the function space 

H{^) := {q G H-^^) : {q, l)_i,i,^ = 0}. 

We also recall that the single layer potential satisfies the following classical jump conditions 


[^^gk = 0, [dn{y<^q)]<^ = F (2.2) 

Let us focus now on the trace of the single layer potential. 

Definition 2.2 For every q G F [~2 (‘^)j we denote by S^q the trace of the single layer operator 5^<^q 
on qf. 


The operator is an integral operator with weakly singular kernel as well. For q G L°°{^) and for 
every x G it reads: 

S'^g(x) = / G(x - y)q{y) dcr^. 


The trace of the single layer operator defines a bounded linear operator from H 2 ('^) into H 2 (^). 
Furthermore, using Green’s formula and the asymptotics (2.1), we can easily prove the identity 


{q,Srq)_i i.g = [ |V(^^g)|^dx <-hoo, Vg G #('^). 
2 ’ 2 ’ Jr 2 


(2.3) 


Following m Theorem 8.15], we also introduce the following particular density and constant which 
will play a crucial role in our analysis. 
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Definition 2.3 The equilibrium density for ^ is the unique density £ H 2 such that is 

constant on ^ and 

= 1 - 

The logarithmic capacity Cap('^) 0/^ is defined as being the positive constant: 

Cap('^) = exp (-27rS<^e<^). 


Setting 

H{^) := {q G hI{^) : (e^,g)_i = 0} 

we know, following McLean m, that the linear operator; 

S^:qG H{^) ^ q G H{^), 

defines an isomorphism that extends into an isomorphism from onto if Cap(^) 1 

(see Theorem 8.16]). Under this condition, one can identify via this isomorphism any density 
q G with the trace 

q := G 

This identification will be systematically used throughout the paper, using the notation with (respec¬ 
tively without) a hat on the density like quantities (respectively on traces). This isomorphism turns 
out to be an isometry provided the spaces H{'^) and Hfif) are endowed with the following inner 
products: 

Definition 2.4 For all q,p ^ Hfia), we set: 


= {tp)-i.i. 


According to (2.3), the inner products introduced in Definition 2.4 are related to the Dirichlet 


energy of the single layer potential through the following identities: 


Uf_lcg = h\\\^= ! |V(^^g)pdx, \/qeH{^). 

We also need in the sequel the following orthogonal projections. 

Definition 2.5 Letll<^ andll<^ denote respectively the orthogonal projections from into H{'W) 

and from into 

In particular, we have the following unique decompositions: 

VgGiL“5(^); q= {q,l)hg+ qQ, %:= U<ffq e H{^), 

VgGiL5(^); q = {e‘^,q)l + qo, qo\= Il<^q £ 


Definition 2.6 We denote by Tr>^ the classical trace operator (valued into H 2 (fif)), and by Tr^ when 
it is left-composed with the orthogonal projection onto HfitF): Tr^ := II-g^Tr^. 
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Let us conclude this preliminary section by a useful characterization of the chosen norm on H{^). 
Classically, we define the quotient weighted Sobolev space; 

VLo^(M 2 ) = {ue P'(M 2 ) : pu E Vu G /R, 

where the weight is given by 

p{x) := ^\/l + |x| 2 log (2 + |xp)^ , X G 

and where the quotient means that functions of LLq (M^) are defined up to an additive constant. This 
space is a Hilbert space once equipped with the inner product: 


Vu-Vudx. 

° Jr2 

In particular, according to (2.3), G ITq^(M^) if and only if g G H(^), and moreover 


Lemma 2.7 For every q G H{^), we have 

= inf {ll^^llwoHR2) : u G and Tr^u = g| . 

The infimum is a minimum which is uniquely achieved by u = d^^q. 

Proof : Given q G let us consider the orthogonal decomposition: 

W^{R^) = 

Let u G lyQ^(M^) be such that Tr^u = q. Writing u in the form u = XS^^q + v with A G M and 
V G {y‘^q)-^, and taking the projected trace on ^ we get: 

q = Xq + Jr^v in 

Forming now the duality product with q and taking into account that: (g, Tr^u)_ i i <^ = {^■^q, ^^)vi/o©k 2 ) 
0, we deduce that A = 1. Since we have now 

= Ikll + IkllwobR^)’ 

with V such that Tr^u = 0, the conclusion follows. □ 
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2.2 Boundary interaction and single layer potential 

In this section, we are interested in quantifying the Dirichlet energy variation between and 
where p = (he. p is the trace of the single layer potential on 7 ). 

Definition 2.8 We define the boundary interaction operators Kp and between T and 7 by: 

\q-.qe H{T) ^ TrO(^rg) e H{^), : p G H{^) ^ Tr° G H{T), 

where Tr° and Trp are given in Definition 

The next result shows that Trp can be replaced by Trr in the definition of 
Lemma 2.9 If p G H{'y), then q := Jrr{<9^-yp) belongs to H{T). 


2.6 


Proof : Let p and q be given as in the statement of the lemma and let us define the function w := 
{wi,We) in by setting: Wi = in hi and We = ^'yP in \ tl. According to (2.3), the 

function w has finite Dirichlet energy since p G H{'y). Thanks to (13), we see that We{x) = 0{\x\ 
at inhnity, and this allows us to obtain the following classical integral representation formula for every 
X G \ hi (see for instance [iQl p. 182] or [TTl Lemma 3.5]) 


We{x) = {dnG{x - ■),We)l r - {dnWe,G{x - ■))_! 1 p 

2 ’ 2 ’ 2 ’ 

0 = {dnG{x - ■),Wi)_i Ip- {dnWi, G{x - •))_ 1 1 p. 

2 ’ 2 ’ 2 ’ 2 ’ 

Since TrrWi = TrrWe = q, we get by subtracting these identities that: 

We{x) = (f, G{x — ■))_i_ 1 r = yTx{x), X G \ hi. 


where the density r := dnWi — dnWe belongs to H(T) since, as already mentioned, w has finite Dirichlet 
energy. Taking the trace on T, we deduce from the above relation that q = r € H{T) and the proof is 
complete. □ 


Proposition 2.10 The operators Kp and K]^ are compact, one-to-one and dense-range operators. 
Moreover, for every functions q G H{T) and p G H{'y), we have: 

= {q, K^p)i p. 


Proof : The compactness follows from the regularity of the single layer potential away from the 
boundary, combined to m Proposition 13.5.8]. 

Addressing the symmetry property, consider q G L°°{T) Ci H(T) and p G L°°{'y) n H{'y). We can 
write that: 



y)q{y)p{x) daydao; 



y)q{x)p{y) dayda^ 


i^^P, q) 


T’ 


and the conclusion follows by density. 
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Assume now that q G H{r) is such that Kp (7 = 0. By the unique continuation property for 
harmonic functions, it means that is constant in and hence that q = 0. Since 

Ran Kr = Ker Kj;, 

we get the density result and the proof is completed. □ 


Proposition 2.11 The norms of the operators Kp and are strictly less that 1. 


Proof : According to Lemma 2.7 we have, for every q G H{T)- 


IKp^lli,.^ = inf |||n||^i(]R 2 ) : u G and Tr°n = Tr°(=yrg)} • 


We deduce that ||Kpg|| i ^ ^ ~ r norm of Kp is no greater than 1. 

The operator Kp being compact, its norm is achieved by some qr G H{T). If ||Kpgr|| i ^ = lkr|| i p) 
we would have, according to Lemma |2.7[ 

-^7^7 = in M^, 


where q^ := Kp^r- This identity implies that gr = [dn{^-yQ'r)]r = 0, yielding the expected contradic¬ 
tion. □ 


2.3 Integral formulation and well-posedness 

Let us go back to Problem Without loss of generality, let us assume from now on that the 

diameter of 11 is less than 1 (otherwise, it suffices to rescale the problem), which implies in particular 
that Cap(r) < 1 and Cap( 7 ) < 1 (see [IHl p. 143] and references therein). 


Proposition 2.12 For every f G HfT), denote by {u^,c^) G H^iPl \ ca) x 
of System (1.2). The function 
follows: 

= F^^p + 


the unique solution 
can be represented as a superposition of single layer potentials as 


w 


(2.4) 


where p G 77(7) o,nd q € H 2 (r) solve the following system of coupled integral equations on the 
boundaries 7 and T.- 


p+ Tr^(yfq) = c-^ on j (2.5a) 

Trri^'yP) + q = f on r. (2.5b) 


Proof : It is a consequence of 133 Theorem 8.16] that the unique solution to System ( |1.2[ ) can be 
written as a superposition of two single layer potentials respectively supported on T and 7 and respec¬ 
tively associated with the densities {p,q) G 77 ~2(7) x as in (2.4). It only remains to verify 

that p is in fact in 77(7), i.e. that {p, l)_i 1 = 0. This is a straightforward consequence of (1.2d) 

and the jump relation for the normal derivative of the single layer potential. □ 








3 The reconstruction formula 


Going back to the DtN operator of problem B, and due to (|1.2d[), we have by Green’s formula 


{dnu^, = -{dnu^, = 0 , 


.7 


which shows that is valued in H{T). Considering data / G H(r), we can thus define the DtN 
operator as follows: 

A^ : / E H{T) ^ dnu^ E H{T). (3.1) 

In the case where a; = 0 , we will denote respectively by Uq and Aq the solution and the DtN A^. 
Note that we have in particular 

Uo = ^r/- 


3.1 Factorization of the DtN map 

Theorem 3.1 The two following bounded linear operators in H(T): 

R := Sr(A.^ - Aq) and K := 
satisfy the following equivalent identities: 

R = (Id-K)-iK, K = (Id + R)~^R. 


(3.2) 


Proof : Giv en / in H{T), let {p,q) E H{'y) x H 2 (r) be the solution of System (2.5). According to 
Lemma 


2.9 


Trr(A^ 7 P) E H{T) and hence Trr(A^ 7 p) = Since / E H{T), we deduce from (2.5b) 

that q = f — y(!yP E H{T)- Applying the projector to (2.5a), we obtain the following system: 


p + Kj.q = 0 
Kf^p + q = f 


on 7 

on r. 


(3.3a) 

(3.3b) 


Eliminating p, it follows that (Id — K)q = / and hence (Id — K)(q — /) = K/. The operator K being a 


contraction (see Propositon 2.11), we end up with: 

g-/ = (Id-K)-i 


K/. 


(3.4) 


On the other hand, we have 

(A 7 - Ao)/ = dni^vq + A^ 7 P)|r - <9n(.5^r/)|r- (3.5) 

But outside P, the two single layer potentials A^r/ and both solve the well-posed Dirichlet 

exterior boundary value problem: 


—Art = 0 


u = f 

u = 0 (|x|“^) 


in \ D, 
on P, 

Ixl —)• -|-oo. 
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Hence {y'rQ + ’^■yP) = =^r/ in \ Q, and in particular we can rewrite (3.5) as 

(A^ - Ao)/ = [dn{yrq + ’^yP)]r - [<9n(^r/)]r = q- f, 


where the last equality follows from the jump relation (2.2). Comparing this relation and (3.4), we 
obtain that Sr(A.^ — Aq)/ = (Id — K)“^K/, which is exactly the first relation in (3.2). The second 
relation follows easily. □ 


The first equation in (3.2) can be seen as a factorization of the (known) DtN operator A.y — Aq 
in terms of the (unknown) boundary interaction operator Kp and Similarly, the second equation 
in (3.2) can be seen as a factorization of the boundary interaction operator K = in terms of 


the measurement operator R (which is entirely determined by the perturbed and unperturbed DtN 


maps and by the exterior boundary T). Using Proposition 2.10, it is worth reformulating this second 
equation in a variational form: 


(K?/, 1 = ((Id + R)-'R/,5) 1 r, V/,5 G H{r). 

i- i- 2 w 2 ’ 


(3.6) 


2 ’ ^ 2 ’ 

This identity constitutes the first step towards the reconstruction of the unknown boundary 7 . Indeed, 
the bilinear form (Kp •, Kp •) 1 ^ turns out to encode the geometry of the inclusion, as shown in the 
next section. 


3.2 Harmonic polynomials and GPST 

Throughout the paper, we identify x = (xi,X 2 ) in with the complex number z = xi -h 1 x 2 - 
Definition 3.2 For every m ^ 1, we define the harmonic polynomials of degree m: 

Pf^ix) = Re (z™) and P^{x) = Im {z^). 

We define as well 

QT,r{x) ■■= Prix) + cT,r and Q^r(^) := PPix) + c^r (3-7) 

where the constant c^p G M, = 1,2, are chosen such that the trace of on T belongs to H{T). 
Finally, we set 


Q? ■■= QtT + iQ2,V 


(3.8) 


The crucial point about these polynomials ^ = 1,2, lies in the fact that since they are harmonic, 
we have 

(Knm I \=Q^, 




,rir^ 




and hence, using these harmonic polynomials Q^p in formula (3.6) (and using for simplicity the same 
notation for the functions and their traces on the boundaries 7 and T), we obtain that for all m,m' 1 
and all i, = 1 , 2 : 


Remark 3.3 The quantities are strongly connected with the so-called Generalized 

' ' u ,'7' 2 iR 

Polya-Szego Tensors (GPST) appearing in the high-order asymptotic expansion of the DtN map for 
small inclusions (see for instance the recent papers by Ammari et al. and references therein). Our 

definition is somehow different from theirs, as they use real polynomials x"^, while we use harmonic 
polynomials. 
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3.3 From the GPST to the geometry of the cavity: an explicit inversion formula 


In this section, we are going to see that the quantities {QT^iQ^k ^ for m, m' ^ 1 and / = 1,2, 

__ 2 ’' 

which can be deduced from the measurements (see (3.6)), contain all the necessary information to 
reconstruct the cavity. We can even say more: the geometric information of 7 is actually redundant 
in the GPST. As we shall see, the knowledge of the quantities and suffices 

to reconstruct the cavity. More precisely, assume that the geometry of 7 is described through the 
conformal mapping 


(j) : z aiz + ao + 


—1 


(Xm^ 


that maps the exterior of the unit disk D onto the exterior of uj (see the book of Pommerenke |38L p. 5] 
for the existence of such a mapping). In particular, t g] — vr, tt] 1 —i^(e**) provides a parameterization 
of 7 . Notice that in this description, |ai| is the logarithmic capacity of 7 and can be chosen such that 
ai > 0. The coefficient oq is the conformal center of uj. With these notation, we have the following 
result. 


Theorem 3.4 Let (Q!^)m^i be the complex harmonic polynomials defined by (3.7)-(3.8). Define the 
two following sequences of complex numbers (1 ^ + 00 ): 


:= = ^(QF,(Id + R)-iRQf)i^r, 

with R := Sr(A.y — Aq). Then, > 0 and we have the explicit formulae: 


oi = 


fii 

27r 


ao — 


T2 

2/ri 




ao 




m 1 , 


where 


and 


Am := {a G : q-q + 2 ai + 802 + ... + (m + l)am = (m + 1), ao m + 1 } 

(_l)|a| + l (27r)T-foi+-+"™) 


: = 


2“om 




(3.9a) 

(3.9b) 

(3.9c) 

(3.9d) 

(3.9e) 

(3.9f) 


The rest of this section is devoted to the proof of this result. 

To simplify the forthcoming computation, we complete the sequence of complex numbers (a^)*,^! 
by setting = 0 for A: ^ 2. We denote a^ (n G N, A: G Z) the A:th coefficients of the Laurent’s series 
of (jf". 

” (3.10) 


— / (Xai^a2 • • ■ ^ani 


\a\=k 


where the sum ranges over all the multi-indices a = (ai, ..., an) G TT whose length |q!| = ai + .. . + ar^ 
is equal to k. We also introduce the quantities: 
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so that 


^m,m' ^ ^ (Q™, (M + R)-l RQ™') , 


IJ-m = and = V 


m,l 


Lemma 3.5 Denoting, for every 1: 


iz) = ^aTz^ and ^ (z) = 


„Tn k 

Ok ^ j 


k>l 


fcsS-1 


the following identities hold true: 

^m,m' ^ r eit( 0 m')/(e*i)<^™(e**)dt and = / e^^cf^')'{e^^)f^{e^^) dt. 

J —TT J —77 


(3.11) 


Proof : Let m,m' ^ 1 and ^ = 1,2 be fixed. For the sake of simplicity, we drop in this proof the 
dependence with respect to 7 and we denote simply by Q™). We aim to compute the quantity: 

(Qr,Qr')i, = (Qr,Qr')-ii,- 


To do so, we recall that from the jump relation (2.2), we have Q™ = [dnU'f^]^, where [7™ := . 

Let us denote by and the restrictions of respectively to \ a; and w. 

We know that solves the following exterior Dirichlet boundary problem: 


-At/-, = 0 


m \ a; 
t/-, = Qr on 7 , 

= 0{\x\~^) as |x| —)■ + 00 . 

The functions u—^ := U^^{4>) are harmonic in \ H (D denotes the unit disk) and satisfy: 

u-i(x) = Re((?i-(^)) + cr and u^ 2 {x) = Im {cf^{z)) + . 

We can easily compute the constants c— and c™ by writing that 

c- + ic^ = - / (Pf^ix) + i P^{x))e^{x) da. = - T <(--(e*')|(/.'(e'*)|e^(e'‘) dt. 


(3.12a) 

(3.12b) 

(3.12c) 


But we know (from direct computations or from [251 Theorem 17.3.3]) that e^{e'^^) = l/(27r|(/)'(e**)|). 
It follows that c— + tc— = —and we have, on the boundary of D: 

<i(^) = ^ [4>^iz) + - Re (a^) and uf^^x) = [^(^) - - Im (a^). 


2 L" 


This can be rewritten, using the identity z = 1/z on dD, as: 


1 




uT, 2 ix) = “2 ) + 4>-{z) - (-(()-(z-i) +(()-(z 


(3.13a) 

(3.13b) 
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These expressions lead us to introduce the following functions: 

W^{z) = + (tr{z) = (Tiz) -a^ + XT(z) (3.14a) 

w^{z) = -^(z-^) + (tr{z) = (Tiz) - a(r + (3.14b) 

where 

\r{z)=W^{z-^)-cl^^{z) and A^(z) =-^(z"!) - (/)-(z). (3.15) 

The functions and are holomorphic in C \ T> and: 

= Re (rc)”) and = I™ (^D 

For every X = (Xi, X2) G identified with Z = Xi + 1 X 2 G C, we have: 

VC/™i(</-(z)) • X = Re [{wT)'{z)Z/ct^'{z)\ and VU^^^ic^iz)) • X = Im [«)'(z)Z/,/.'(z)] . 

On 7, the outer unit normal vector is parameterized by t G [—7r,7r[i—>■ e*V^(e**)/|</>^(e**)|, and therefore, 
for every m' ^ 1: 


/ dnU:^,{x)Qf{x)da,= / Re [e''«)'(e'*)] Re {e^^)]dt. 
J'y J — TT 

On the other hand, solves the following interior problem: 


-Aun = 0 


m ui 


U^, = QT on 7, 

whose unique solution is merely C/j™ = QJ^, so that: 

/ dnUr;[{x)Qf{x)da,= r Re [e*‘(,/.”^)'(e*‘)] Re [r'(e**)]dt. 


Gathering now (3.16), (3.17) and taking into account the expressions (3.14), we infer that: 


{QT,QT')l^ = J [dnUr{x)]^QT'{x)da, = - j Re [e'*(Ar)'(e'‘)] Re [</.™'(e'*)]dt. 
Notice now that, on dD, we have: 

Re [e**(A)")'(e**)] = - [e**(<^™)'(e**) + e*'((/)™)'(e*')j = -2Re [e**(</>((^)'(e**)] , 
and therefore, (3.18|) can be rewritten as: 


(Qr,Qr)i^. = 2 




Re [e**(<^™)'(e**)] Re U™'(e**) 


Using similar arguments, lengthy but straightforward computations lead to: 


{Q^,Qf)i =2 I Re[e^\cj>^y{e^^)]lm 0™'(e'*) 


= 2 / Im [e**((/.™)'(e**)] Re 0™'(e**) 


dt. 

dt 

dt 


(3.16) 


(3.17) 


(3.18) 
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and 


= 2 r Im [e**(<^™)'(e**)] Im U-'(e**) 
J —TT ^ 


dt. 


Formulae (3.11) follow. 


□ 


Lemma 3.6 The following relations hold true: 

, T2 

0-1 = \ — and ao = -—. 

V 27r 2/xi 

Proof : For every z ^ C\D (recall that D denotes the unit disk), we have <fi\{z) = (j)+{z) = aiz and 
4>\{z) = {ai)‘^z‘^ + 2aiao^;, and hence 

z(l)'_^_{z) = aiz and z{(f)‘^y{z) = 2{ai)‘^z‘^ + 2aiOo. 


Applying formulae (3.11), we obtain: 

fii = 27r(ai)^ and fj .2 = 47r(ai)^ao. 

The conclusion of the lemma follows. 

The conformal mapping (j)~^ can be expanded as a Laurent’s series having the form; 

biz + feo + ^ hz^, 


□ 


/.-i, 


z = 


outside a disk D' centered at the origin and containing uj. The complex coefficients bk {k ^ 1) can be 
deduced on the one hand from the coefficients of (/>, and on the other hand from the values of t'm, 
{m ^ 1), as claimed in the following lemma: 


Lemma 3.7 The following relations hold true: 

bi = 1/ai and 6o = —ao/ai. 

For every m ^ 1, we have: 


b—m, — 


27raim m 


- E 

m < ^ 


ai3i ■ ■ ■ ap ^, m ^ 1. 


(3.19) 


(3.20) 


1=-i 


Proof : Identities (3.19) follow straightforwardly because (f and (j) ^ are inverse mappings. Integrating 
by part the expression of Um = obtained from (3.11), we get: 


n^ = ai / e*V”'(e**) dt = -aim / e*V'(e*‘)e*V”'" (e**) dt 


(3.21) 
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Since 1 1 —>■ is a parameterization of dD, applying Cauchy’s integral formula we get on the one hand 


/ e*V"*(e**) dt = -i / = 2™™!- 

'—TT J dD 


On the other hand, since 1 is a parameterization of 7 and the function (p ^ being holomorphic 
in C \ cj, we have: 


J —TT 


= -i / de = 27r6_, 


IdD' 


Identity (3.21) can thus be rewritten as: 

Vm = 27raia™i = —2TTaimb-m, 


and identity (3.20) follows according to (3.10). 


□ 


Using the above lemmas, we are in position to prove the main result of this section, namely 
Theorem 13.41 


Proof of Theorem 3-4' Since cj) and (f ^ play symmetric roles, we can exchange am and bm in Formula 
(3.20) to obtain: 


m — 

m 


I E 


'/3i • • ■ 


m ^ 1. 


1=-i 


Reordering the terms of the above sum, we get that 


O—m — 

m 


1 E 


(3.22) 


(0,a)eB™ 

where is the set of {0, a) G N x such that 


0 + ao + «i H-1- ccm = "i 

6 - {ai + 2a2 H-h mam) = -1- 


Now, one can easily check that {0,a) G if and only if a belongs to the set Am defined by (3.9e) 
and 6 = m — (ao + ai + • • • + am)- Therefore, (3.22) also reads 


^—m 


1 

m 




Using (3.19) and the hrst equality of (|3.20[) in the above relation, we obtain that 


^—m 


1 

- a 

m 


—m 

1 


E(-»or(-|^)“‘- 

CuSAm 


V 27rm) ’ 


and the conclusion follows immediately. 


□ 
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3.4 About stability 


It is well-known that logarithmic stability is best possible for Calderon’s inverse problem. In the 
particular case of cavities, this result is proved in [21 Theorem 4.1] where the error on the geometry 
(measured using the Hausdorff distance) is estimated in terms of the error of the DtN (measured in 
operator norm). However, as suggested by Alessandrini and Vessella [3j, one can try to construct 
stable functionals, namely Lipschitz-continuous functions of the data carrying relevant information on 
the geometry of the obstacle. According to formula (3.9), each coefficient a^. A: ^ 1, yields an example 
of such functional. Actually, we can prove that each coefficient is not only a Lipschitz-continuous 
function of the data, but is analytic. Let us define the following open subspace of C{H{T)): 


= {R £ A2(LI(r)) : Id -I- R invertible and /ri(R) > 0}, 


where /Ui(R) := (Qp, (Id -|- R)~^R(5p) i p. 

Notice in particular that, for every Lipschitz Jordan curve 7 , the continuous linear mapping R = 
Sr(A.y — Aq) belongs to Up- We deduce straightforwardly the following analyticity result: 


Theorem 3.8 On the open subset Up of C{H{T)) define the sequence of analytic functions ak'-Up ^ 
C (k ^ 1) as given by the formulae (3.9). // R = Sr(A.y — Aq) for some Lipschitz Jordan curve 7 , a 
parameterization of 7 is given by: 


t e] — vr, tt] ^ Ofc(R)e*^* G C. 
fcsgl 


4 Numerical results 


We present in this section some numerical experiments meant to illustrate the feasibility of the pro¬ 
posed reconstruction method. For the sake of clarity, let us first sum up the different steps of the 
simple reconstruction algorithm: 


1. Compute a numerical approximation of the operator R = Sr(A-^ — Aq). 

2. Fix an integer M and compute for 1 ^ m ^ M 


hm ■■= = ^(gF,(Id + R)-'RQf)i,r 

:= = ^(QF,(Id + R)-'RQf)i,r. 


3. 

4. 


Compute (a_m)-isgm 7 M via formulae (3.9). 
Plot the image of the unit circle by 


(pMiz) = aiz + ao + 


E 




—m 


l^m^M 


Let us give some details about the implementation. We use the finite dimensional approximation space 
spanned by the family := {Q(?, Q™, 1 ^ m ^ M}. We denote by Qr the 2Mx2M complex matrix 
whose entries are the {f,g)i p, where (/, 5 ) G Qp x . Note that Qr is nothing but the Generalized 
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Polya-Szego Tensor (GPST) associated to P. Obviously, a similar matrix can be defined for the 
boundary 7. We denote by R the matrix whose entries are (/, Rs')i p = ((A-v — Ao) 5 ,/)_i 1 p, for 


(/)<?) £ Qr ^ 2r ■ With this notation, the reader can easily check that formula (3.6) admits the 
following discrete version: 


Q.y ~ Q.y := Qr(Qr + R) ^R- 


(4.1) 


This formula relates in a very simple way, through the measurement operator R, the GPST of 7 to 
the GPST of P. In particular, the coefficients Hm and are particular entries of 5 Q 7 . 

We consider now a test configuration in which P is an ellipse centered at the origin and of major 
axis [—1.9,1.9] and minor axis [—1.1,1.1]. The boundary 7 of the obstacle is parameterized by: 


1 

t G] - 7 r, 7 r] ^ aA;e*^^ 

k=-7 


where the complex coefficients are given in the following table: 


ai 

Oo 

0-1 

0-2 

0-3 

a_4 

0-5 

0-6 

0-7 

0.5 

-1 

0.085 

-0.06i 

-0.035 

0.06i 

0 

-O.Olf 

-0.005 


The data are generated using the Matlab Laplace boundary integral equation solver (for more 
information, see this link: lES). Taking M = 12, we hrst show on Figure]^ the reconstructed cavity 
for exact data and using the eight coefficients oi,..., a_ 6 . 



Figure 2: Typical configuration: reconstruction with ai,...,a _6 (in red) and actual inclusion (in 
gray). The blue point stands for the position of the origin r. 


Instead of using the harmonic polynomials z” in Definition 3.2, one can use the shifted harmonic 
polynomials (z —r)”, for some given r G C. This additional parameter turns out to have some influence 
on the quality of the reconstructed cavities, as shown in Figure [^ For instance, choosing r in the 
neighborhood of —0.5, one can recover the six coefficients ai,... ,a _4 with a relative error less than 
2 %, while for r = 0 this accuracy is achieved only for the coefficients ai,..., a_ 2 . 
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Figure 3: Relative error of the retrieved coefficients (in %) with respect to the abscissa of r (there is 
no relative error for the coefficient o _5 because it is null). 


Let us take now r = —0.5 and consider a more realistic configuration of noisy data. We generate a 
random matrix N having the same size as R and whose coefficients are uniformly distributed between 
—1 and 1. For 5 = 0.05,0.15,0.25,0.35, we compute the matrix whose coefficients are: 

Rfj = (1 + 5Nij)Rij , 1 ^ i, j ^ 2M, 

and we replace R by R^ in formula ( |4.1[ ). 

We show on figures |4]j^ examples of reconstructed cavities respectively with 5%, 15%, 25% and 35% 
of noise. The number of correctly recovered coefficients decreases with the level of noise and only those 
coefficients are used in the reconstruction. We plot on Figure the dependence of the mean relative 
error with respect to the level of noise and we notice a good stability of the first three coefficients 
0 . 1 , oq and a_i. 

Finally, we illustrate on Figure]^ the efficiency of the method for more complex geometries (non 
convex outer boundary F and a non centered cavity). The choice of the parameter r to obtain good 
reconstructions is not clear so far and this would need to be further investigated. 


A Appendix 

Consider problem 0 with a piecewise conductivity (t{x) = 1 + (a — l)lu]{x) (where l^j denotes the 
characteristic function of oj and a a positive constant). 


Proposition A.l For every f G H 2 (r), System (1.2) admits a unique solution (tt-^,c-^) G H^{Q\ 
w) X M. It is the unique pair that realizes: 


mm 

{u,c)&m{n\Lj)xi 


/ |Vttpdx : u|r = /, u\^ = c 

J Q\cJ 


(A.l) 


The function can be considered as a function of H^{Q) by setting = c^ in u. 
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1.2 



Figure 4: Reconstruction (in red) with oi,... ,a _4 and actual inclusion (in gray) with 5% noise. 


For every a > 0, System (1.1) admits a unique solution uh G This function achieves: 

]- [ \Vu\^dx + ^ f iV-updx : u|r = /l. 

^ J ^ Juj I 


mm 

u&H\n) 1 2 


(A.2) 


The following convergence result holds true for every f G 

—)• u-^ in Ff^(n) as a ^ + 00 . 


Proof : The minimization problem (A.l) can be reformulated as: 

|V(rc + e^ + cv)fdx, 


mm 

{w,c)£H^{n\uj)xRjQ\ij 


(A.3) 


where e-^ and v are both harmonic in n\a; with Dirichlet data e-^|r = /, = 0 and u|r = 0, v\^ = 1. 

For every w G \ ^)i we get: 

/ \'V{w + C'^ + cv)\‘^dx = / |Vuipdx+ / |Ve'^pdx + c^ / |Vupdx + 2c / Sne-^du, 

jQ\iJ JD.\cJ JD.\uJ J'y 


and therefore the minimum in (A.3) is unique and achieved for w = 0 and 

c^ = — f dne^dcri f |Vu|^dx) 

J'y \^Jq\(U J 

The corresponding fonction u^ ■.= e^ + c^v is the unique minimizer of problem (A.l), and can easily 


be shown to solve System (1.2). It is classical to verify that, reciprocally, every solution of System 


(1.2) provides a solution to the minimization problem (A.l) 
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Figure 5: Reconstruction (in red) with oi,..., a _4 and actual inclusion (in gray) with 15% noise. 


Seeking the minimum of problem (A.2) in the form u = + w with w G we are led to 

consider the new, equivalent, minimization problem: 


mm 

w£H, 


in / iVrcpda; + ^ / |Vt(;pdx+ / dnU'^'wda, 
^ Jto J'y 


(A.4) 


where we have used the fact that in a;. The existence and uniqueness for such a problem is 

straightforward and we denote by Wa the minimizer. Introducing 


—f 
W'^ = 


1 


mes(a;) 


w. 


idx, 


and taking into account the condition (1.2d), the last term in the right hand side can be rewritten as: 


/ dnU-^w^da= / dnU^{wl^ —wl) da. 
f "y J 'y 


Invoking the Poincare-Wirtinger inequality in uj, we get the estimate: 


dnU^wl^ dcr 


(A.5) 


where the constant C > 0 depends only upon oj. The minimum (A.4) is negative (tc = 0 is an 
admissible function), whence we deduce that: 


<C, 


and therefore 


0 as a-^+oo. (A.6) 

Remarking again that the minimum (|A.4[) is negative and using the estimate (A.5) together with the 


convergence resnlt (|A.6|), we deduce that: 

as a —+ 00 , 
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Figure 6: Reconstruction (in red) with oi,..., a _2 and actual inclusion (in gray) with 25% noise. 


and the proof is completed. 


□ 
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